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ABSTRACT 

We carried out an unbiased, spectroscopic survey using the low-resolution module of the infrared 
spectrograph (IRS) on board Spitzer targeting two 2.6 square arcminute regions in the GOODS-North 
field. IRS was used in spectral mapping mode with 5 hours of effective integration time per pixel. 
One region was covered between 14 and 21/zm and the other between 20 and 35/irn. We extracted 
spectra for 45 sources. About 84% of the sources have reported detections by GOODS at 24/im, with 
a median f^{2A^m) ^ IQQjiJy. All but one source are detected in all four IRAC bands, 3.6 to 8 
^n\. We use a new cross-correlation technique to measure redshifts and estimate IRS spectral types; 
this was successful for ~ 60% of the spectra. Fourteen sources show significant PAH emission, four 
mostly SiO absorption, eight present mixed spectral signatures (low PAH and/or SiO) and two show 
a single line in emission. For the remaining 17, no spectral features were detected. Redshifts range 
from z ~ 0.2 to z ~ 2.2, with a median of 1. IR Luminosities are roughly estimated from 24/im flux 
densities, and have median values of 2.2 x lO^^Lo and 7.5 x lO^^L© at z ^ 1 and z ^ 2 respectively. 
This sample has fewer AGN than previous faint samples observed with IRS, which we attribute to 
the fainter luminosities reached here. 
Subject headings: galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 

The Spitzer Space Telescope (Werner et al."2004) has 
multiplied by several orders of magnitude the volumes 
previ ously surveyed in th e infrared for extragalactic ob- 



jects (Soifer et al. 20081. Mid-infrared continuum sur 
veys turned out to be especially fertile ground because 
of the superb sensitivity and speed of the 24/i m chan- 
nel in the MIPS instrument ( |Rieke et aL]|2004[ ). These 
mid-infrared surveys have established the dramatic evo- 
lution of g alaxy populations over the la st 10-12 billion 



years (e.g. |Le Fioc'h et al.| |2004[ |2005| ; [Caputi et al 



et al. 2001 



2006 PapoviclT'eFar^OT') previously suggested by ISO 
surveys~(e.g . |Elbaz et al i""" '^i- ^ r,nnn t^^i_ 



^ 1999: Ohver et al.. 2000, , Dole, 

^ ^ ,owan-Robinson et al. 2004|), and yielded 

rich samples oF galaxies for followup. Critical to the in- 
terpretation of the data is the availability of multiple- 
band detections, which are used to constrain the red- 
shifts, luminosities and sometimes the powering mecha- 
nism of sources. These constraints come about because 
the mid-inf rared spectra of galaxies are often str ongly 
structured ( Smith et al. 2007b Armus et al. 2007 ) , and 
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and to heating and dust geometry. By the same token, 
any single-band survey will have redshift-dependent bi- 
ases resulting in a sampling of the population that varies 
with redshift. The obvious way to avoid this sampling 
bias is to survey in the dispersed light, a technique pi- 
oneered from the ground using objective-prism imaging 
in the visible. 

In this paper we report on a spectrally unbiased extra- 
galacti c census in the mid- infrared using the IRS instru- 
ment (Houck et al. 2004) on Spitzer, and covering two 
spectral windows, 14 to 21/im and 20 to 38/im. Besides 
the specific results reported here, this survey illustrates 
the advantages of spectral mapping in the infrared. The 
placement of the windows allows us to attach an observed 
spectrum to the corresponding 24/im source, and could 
therefore potentially address the bolometric correction 
for each source. The data also allow us to use spectral 
shapes, features and emission lines to characterize all de- 
tected sources, rather than only those sources selected 
for follow-up by the inevitably rough criteria of filtering 
on colors within the mid-infrared or in combination with 



other bands (e.g. Yan et al. 2004 
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examination of the emission line properties also yields a 
fair assessment of the distribution of line strengths and 
line-to-continuum ratios. Such an assessment constrains 
the frequency of line-dominated sources, and might even 
yield examples of sources radiating largely line emission, 
and therefore very rarely picked up in broad-band sur- 
veys. When our survey was proposed, a few examples of 
line-dominated source s in the m id-infrared were known, 
including NGC 1569 ( |Lu et al.| |2003). Today the most 
striking case known is probably the intergalactic shock 
in Stephan's Quintet, which emits more than 20% of its 
total infrared emission in the main three molecular hy- 



drogen pure rotational transitions ( [Appleton et aL]2006 1. 
There is no reason to rule out scaled up versions of these 
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systems, which might be detectable with Spitzer. 

Source confusion is often a hmiting factor for infrared 
surveys, since the need to cool telescopes limits the size of 
the primary. The confusion limit appears in a variety of 
forms and definitions (Hclou & Bcichman 1990: Rowan- 
Robinson|[2001: Lagachc ct al. 20U3 , Dole ct al., 2U04} , 
but is ultimately dictated by the source density and its 
dependence on flux. A spectrally dispersed survey offers 
the spectral dimension as a discriminator among closely 
spaced sources, since the strongly featured shapes of mid- 
infrared spectra would cause sources to peak almost al- 
ways at different wavelengths. A method based on PSF 
fitting at different wavelengths allows to estimate accu- 
rately the confused sources positions at their peak wave- 
length as well as their contribution in flux to the neigh- 
boring sources. This spectral separation, solely relying 
on the knowledge of the PSF of the instrument, can be 
used as a prior to extract distinct spectra and thus disen- 
tangle sources, even if their broad-band images overlap 
substantially. 

This paper describes the survey design, observations 
and data reduction in Sect. 2, describes the extraction 
of redshifts and spectral types from the data in Sect. 
3, presents results on the distribution of redshifts and 
relation to broadband surveys in Sect. 4, and discusses 
those results in Sect. 5. 

Throughout this paper, we assume a A-CDM cosmol- 
ogy with Ho = 70 km s~^ Mpc~^, flM = 0.3 and 
r^A = 0.7. 



2. OBSERVATIONS AND DATA REDUCTION 

The goal was to carry out a blind survey rather than 
targeting specific sources. However, the general location 
of the surveyed area on the sky could be optimized to 
minimize sky brightness, and improve access to ancillary 
data. We chose to place the survey in the Hubble Deep 
Field North (HDF-N) covered in the infrared by Spitze r 



as part of the GOODS survey (Dickinson et al.|2003l 



The availability of deep imaging with Spitzer allows us 
to relate the IRS data to continuum imaging immedi- 
ately, and to avoid bright sources that would generate 
artifacts and reduce the effective size of the survey area. 
Additional ancillary data such as redshifts, ground-based 
photometry or high-resolution imaging all facilitate in- 
terpretation. We also made an effort to place the survey 
area in a particularly low-density region of GOODS-N 
judging by the 24/im continuum map (FigfTl). 

The GOODS 24/j,m data are 84% cornplete around 
SOyitJy, and go as deep as 10/iJy, the weakest extractions 
reported by the GOODS team. 

2.1. Observations 

We observed a region in the Hubble Deep Field North 
(Fig [l]) with the Long-Low module of the Infrared Spec- 
trograph (IRS) onboard Spitzer. Data are collected si- 
multaneously from the two slits (LLl and LL2) of the low 
resolution module of the IRS, which cover different spec- 
tral ranges. We chose to maximize the depth of coverage 
with each slit, and therefore obtained two spatially and 
spectrally disjoint sets of data. The survey was executed 
using 8 AORs (Astronomical Observation Requests) for 
a total of 46.5 hours of observation resulting in an inte- 
gration time at full depth of 5 hours per sky position. 
Each AOR consists of 65 slit exposures following a 13 by 



5 raster map with a step of one pixel (5.1") between each 
exposure. Moreover, the 8 AORs are dithered by 1/3 of 
a pixel in order to improve the Point Spread Function 
(PSF) coverage at short wavelength. Those observations 
yield two adjacent areas, covering 14 to 21/im (LL2) and 
20 to 38/im (LLl) respectively. Each field is covered to 
near-full depth over an area of 165" x 56" (innermost 
contour in Fig [T]) , corresponding to a number of spatial 
resolution elements of 33 x 11 for both LLl and LL2. 
This is the most effective use of telescope time and in- 
strument for an exploratory survey, rather than attempt- 
ing to cover the same area with both sub-modules to get 
complete wavelength coverage. 

2.2. Data filtering and Cube construction 

Processing was pushed farther than the usual IRS 
pipeline products to obtain the best sensitivity. We ap- 
plied additional reduction steps starting with the Basic 
Calibrated Data (BCD) level processed through the ver- 
sion S15 of the pipeline. BCD are dispersed slit images 
that have been cleaned of radiation hit artifacts , rec- 
tified, and calibrated (see IRS data handbook: http : 
V/ssc. spitzer . calt ech . edu/ irs/ dh/ dh32 . pdf ) . 

l^'irst, we corrected for an upward drift in the BCD, 
manifested as a monotonic increase in the median sig- 
nal with time for each set of consecutive AORs. This 
effect was assumed to be a detector effect, since the rise 
is too steep for any plausible phenomenon on the sky. It 
was thus removed as a function of time (robust polyno- 
mial fit). The robust fit was used to be resistant to a 
small number of outliers in the BCD offsets. The pro- 
gram used was robust_poly-fit, available as part of the 
IDL Astro Library. This procedure also removes any real 
time-invariant background, since we end up with a zero- 
median dataset. The data at this point contain solely 
sources plus rogue signals. The latter originate in rogue 
pixels, whose dark current is abnormally high and varies 
with time and different sky brightness. Considering the 
impact of such artifacts on the faint sources we want to 
extract, a two-step removal was applied. The first step 
relies on sets of zero-median BCDs taken close in time 
within a AOR, which we call an ensemble. After choosing 
the length of the zero-median ensemble, rogue pixels were 
subtracted by removing pixel-by-pixel a running robust 
trimmed linear fit to the ensemble from each BCD. The 
program used was robustJinefit, available as part of the 
IDL Astro Library. A robust linear fit rather than a ro- 
bust mean was required to take into account some wave- 
length dependence in the drift described above. However 
this residual drift was small and the fits were close to ro- 
bust means. This method only removes rogue signals 
thus ensuring that no correlated-noise is added to the 
data. We used IRSCLEAN ^ after the ensemble subtrac- 
tion to check how many additional bad pixels are found 
for ensembles of decreasing size. Changes in the num- 
ber of bad-pixels detected were found to be less than 
2% when going below 1/4 of an AOR (32 frames) which 
was then chosen as the ensemble size. Additional rogue 
pixels were found and repaired at the single-BCD level 
using IRSCLEAN. IRSCLEAN fixes pixels by interpolat- 

* software provided by t he SSC. 
http:// SSC .spitzer .caltec h.edu / archanaly / contributed / irsclean / | 
IKSULEAiN JVlASK.htmT] 
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Fig. 1. — Depth-of-coverage of the Spiteer Unbiased Ultradeep Spectroscopic Survey overplotted on MIPS 24/ir?i observation (background) 
of the GOODS-North field (Dickinson ct al. 2003). We show depth contours in white at values of 48, 95 and 143 (number of times a given 
position is covered by the 1KB shtj, respectively trom the innermost to the outermost contour. 



ing adjacent rows. One fixed iteration was done on each 
frame at a threshold of 4 cr. 

The second step removes the remaining rogue pixels 
in the map-making. We used CUBISM developed by 



Smith et al. (2007a I, a custom tool created for the as- 
sembly and analysis of spectral cubes from IRS spectral 
maps. Zero-median BCD frames were recombined into 
two spectral cubes providing 2 spatial and 1 spectral di- 
mensions for each of LLl (Field 1) and LL2 (Field 2). 
CUBISM uses trimmed averages over the set of individ- 
ual measurements going into the signal for each pixel in 
the cube. The measurements are first weighted by the 
overlap fraction between BCD pixel area and cube pixel 
area, then averaged with trimming. We represent one of 
the cubes on a diagram in Fig[2j The considerable redun- 
dancy in the spatial coverage and the dithering between 
AORs oversample the sky sufficiently that the data cubes 
could be built with (2.55")^ pixels, oversampling by a fac- 
tor of ^2 the native spatial pixels. The cubes were not 
oversampled in the spectral dimension, thus the sampling 
remains 0.092/im/pixe/ in LL2 and ranges from 0.177 to 
0.187 fim/ pixel in LLl. The resulting cube dimensions 
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Fig. 2. — Representation of a reconstructed datacube using the 
CUBISM software developped by the SINGS team (Smith et al. 
2007). 



are 79 x 28 x 95 pixels covering 201" x 71" x 17.3/im for 
Field 1 and 79 x 28 x 75 pixels covering 201"x71"x6.9^m 
for Field 2, corresponding to the outermost contours in 
Fig[l] 

2.3. Characterization of the Noise 
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Fig. 3. — Standard deviation of the noise in our dataset in each 
wavelength bin and for both fields (Fieldl in red and Field 2 in 
blue). Solid lines illustrate the noise relative to each source extrac- 
tion whereas the dotted lines correspond to the noise in one single 
pixel of (2.55")^. 

In order to estimate the noise in the datacubes and 
extract its variation with wavelength, we integrated each 
cube over its whole spectral range creating 2D map-like 
planes (one per field). We applied an iterative sigma- 
clipping on those maps to mask pixels containing sources, 
leaving about a thousand pixels in each field, excluding 
the edges where the noise increases dramatically due to 
decreasing coverage in the map. These remaining pix- 
els are assumed to be dominated by noise. Spectra at 
all of these pixels were then used to compute the Icr- 
deviation of the noise at each wavelength. Note that 
due to the background substraction, these spectra have 
a zero mean. Integrated over one oversampled (2.55")^ 
pixel, this standard deviation ranges between 2.7/iJy and 
9.1/iJy in Field 1 (LLl, between 20 and 35/im) and be- 
tween 2.3yLtJy and 6.3/L(Jy in Field 2 (LL2, between 14 
and 21/im) (see dashed lines in Figjs]). The samples of 
positions used to extract the noise were then split into 
several spatially distinct subsets to check for any spatial 
dependency but no significant variation was found. 

All information at A > 35/Ltm in Field 1 (LLl) was dis- 
carded due to very high red-end noise. In the following, 
LLl spectral range refers to a 20 to 35/im band. 

2.4. Source detection and extraction 

The three dimensional structure of the data was used 
to detect sources in both spectral and spatial dimensions. 
Sources were selected down to a low significance by scan- 
ning through X-A planes by eye and then noting their 
positions in X-Y planes (summed over a portion of the 
spectral range) as can be seen in Fig |4] The signal was 
then estimated in the spatial vicinity of each position 
and summed over the spectral range where the source 
signal is higher than 2a, using optimal extraction (Nar- 
ron et al. 20071, applying a matched filter to the data, 
i'his estirnated signal was then fitted by a 2-dimensionnal 
Gaussian function in order to get the best positions and 
the best spectra for our sources. Finally, we estimate the 
integrated SNR for each optimally-extracted spectrum 
over the full spectral range of LLl or LL2 and keep only 
sources with integrated SNR greater than 2. We were 
able to extract 45 spectra, 20 in LLl and 25 in LL2 (see 
Table Jl] and Fig |5] and [6]) . We show in Fig [t] the Signal- 
to-Noise ratios at peak achieved in our spectra. We also 
show the contracted bandwidth in /zm over which our 



spectra achieve a SNR greater than 2. It should be noted 
that in estimating signal to noise ratios, the noise values 
presented in the previous section need to be multiplied 
by a factor of 12 to account for the spatial integration 
implicit in optimal extraction (see solid lines in Fig[3|. 

2.5. Comparison with 24^^m GOODS catalog 

We correlated our sam ple with the M IPS 24/im catalog 
from the GOODS survey Dickinson et al. (2003). Source 
detection in the MIPS image uses IRAC positions as pri- 
ors, thus allowing to select sources to a lower level than 
expected in the presence of confusion. The final M IPS 
catalog is described in more detail by Chary ( |2007| ; it 
has a la depth of about 5/j.Jy and reaches a 84% com- 
pleteness limit of 80/iJy. We select the nearest MIPS 
source up to 3 arcseconds (corresponding to ^ la in IRS 
positioning accuracy) to be the counterpart of our IRS 
detections. We observe a mean difference of l.Tbetween 
IRS and MIPS counterpart coordinates. Results of the 
cross-identification (together with the IRAC catalog) are 
presented in Table [2j Considering both LLl and LL2, 38 
out of 45 sources (84%) of our sample have a 24/im coun- 
terpart with a median f^(24:fim) of 103/iJy. For LLl 
only, the match rate is 90% (18 out of 20 sources). All 
MIPS sources with ^24 > 80/i Jy that fall into the region 
of near-full depth coverage of this survey are detected. 
The faintest MIPS counterpart to one of our extractions 
(SUUSS 18) has a 24/im flux density of 47 /iJy. This 
source is detected in LLl with a maximum SNR of ~7 
at 30.5/im. 

The 24/im fluxes were used to check the flux calibration 
of our dataset on sources extracted from Field 1 where 
the spectral range overlaps the MIPS 24/im filter. We in- 
tegrated the IRS LLl spectra under the 24/im bandpass. 
Sources with SNR greater than 2 in this integration were 
compared to their GOODS 24/im flux counterpart in Fig 
[8) A linear regression on the data in Fig |8] yields 

S24JRS = 1.14 (±0.09) X S24,MiPS - 27.1 (±16.2) 

As one might expect, the main source of error in esti- 
mating the 24/im band fluxes using IRS spectra is the 
low SNR. We don't observe a clear bias or offset in the 
calibration. 

In order to further compare our data to the photo- 
metric Mid-IR data provided by the GOODS survey, we 
convolved the whole Field 1 cube by the MIPS 24/im fil- 
ter. Contours were extracted of the resulting IRS 24/im 
map and overplotted on the released MIPS observation 
at the same wavelength (Fig [9]) . No bias or offset was 
observed in the comparison ofsource positions. We can 
detect a few examples of spatial confusion in the IRS 
map that are not present in the MIPS observation due 
to the better spatial resolution. These confused sources 
were however separated spectrally before extraction. 

3. DETERMINING REDSHIFTS AND SPECTRAL TYPES 
FROM IRS SPECTRA 

We have developed an original method to determine 
the redshift and spectral type of our sources using tem- 
plate spectra typical of various types of galaxies, and 
cross-correlation to fit the observed spectra to the tem- 
plates. We describe here the method and present the 
results. 
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Fig. 4. — On the left, spatial map (X-Y plane, cf Fig|2]| of the Field 2 cube integrated over 14-19/xm.On the right, spatial map (X-Y 
plane) of the Field 1 cube integrated over 21-3Afj,m. Both maps are projected onto the sky. 



We selected a set of 21 template spectra: 5 templates 
dominated b y Arom atic Feature emission prese nted in 
Smith et all ( |2007b| ), 13 ULIRGs spectra from |Armus 
et al. "i [2ijU7p , two radio galaxy and quasar spectra (1^. 
(Jgle private communication) and the spectru m of the 
Wolf-Rayet galaxy NGC1569 f Wu et al.|[2006| . For all 
those templates we have full IRS low resoiution spectral 
coverage from 5 to 38 /im. All spectra were converted to 
restframe before being used. They cover a full range of 
source properties from star-forming to AGN-dominated 
galaxies, and include the various known mid-IR sig- 
natures: Aromatic Features, silicate absorption, both 
high and low ionization lines, non-thermal continuum 
emi ssio n, and steeply rising thermal dust emission (see 
FiglToland Table Isf. 



product, ps,t{z), for that value of z as in Eq. IT] above, 
and refer to it as the cross-correlation function. We build 
cross-correlation functions for < 2; < 3 with Az — 0.01. 
Such functions vary between 1 (perfect linear correla- 
tion) and -1 (perfect linear anti-correlation). A cross- 
correlation maximum indicates a redshift identification 
with high likelihood while the cross-correlation function's 
behaviour away from the peak reflects the shape of the 
template in comparison to the observed spectrum. 

3.2. Effect of Noise in the Spectra on Cross-correlation 

The noise in the data has an important effect on the 
cross-correlation. Assuming that the templates are virtu- 
ally noiseless in comparison to our spectra the measured 
Pearson product-moment correlation coefficient is in fact 



3.1. Correlation analysis of individual sources 

The method is based on a 2 parameter (redshift and 
spectral template) cross-correlation to estimate IRS spec- 
troscopic redshift and determine a best fit spectral type 
for our sources. For a given redshift z, we compute the 
Pearson product-moment correlation coefficient between 
data in the observed frame and each template spectrum, 
redshifted to z. The Pearson product is an estimate of 
the degree of linear relationship between two data vec- 
tors, usually noted ps,T, 



Ps,t{z) = Ps,t{z) X 1 -I- 



<3l 



(2) 



PS,Tiz) 



Cov{S, T) 



(1) 



where cr at is the standard deviation of the noise in our 
dataset (see Fig [3]). This equation shows that the noise 
will produce an overall decrease of the amplitude of the 
cross-correlation function. To quantify the effect of the 
noise on the cross-correlation function ps,t{z) we simu- 
late data spectrum using the templates. The noise in our 
data has been characterized in Sect. 2.3 and its standard 
deviation, (Tat, extracted as a function of wavelength for 
both fields (LLl and LL2) is shown in Fig |3] We define 
the SNR in our spectrum as 



SNR^ 



<S> 



AXmax 



where S and T are the data spectrum and the template 
spectrum, and as and ctt are the standard deviations of 
the source spectrum and the template spectrum in the 
wavelength range (LLl or LL2) respectively. Note that 
as is not the standard deviation of the noise presented in 
Sect. 2.3, but an estimate of the departure of each spec- 
trum from a constant flux density (i.e. the usual deflni- 
tion of a standard deviation if we consider S as a random 
variable). Instead of a traditional cross-correlation with 
lag, we redshift the template by z, compute the Pearson 



< >AAri 



(3) 



where AXmax is a small wavelength band around the 
maximum of the spectrum (2p,m and 1/im in LLl 
and LL2 respectively). We add noise in the tem- 
plate spectrum to reach a speciflc SNR then com- 
pute Pt{z),t{z) (Template-with-noise-added vs Tem- 
plate) which is equivalent to an auto-correlation func- 
tion. We then extract the real value of the peak of the 
auto-correlation and plot it against the SNR in the spec- 
trum (see Fig 11 ). At each given SNR we make 30 real- 
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Fig. 5. — Spectra of the 18 sources extracted from FIELD 1. In shaded gray we plot the Itr-deviation for each spectrum. Two additional 
source spectra were extracted from this field. They are discussed in more details in sect. 4.3. 



isations of the noise and average them to get the mean 
auto-correlation peak value. This relation, built on the 
template, is very helpful to estimate the SNR in our data 
from the value of the peak of the cross-correlation. As we 
can see from Eq. [2] crs/crjv controls the dynamic range 
of the cross-correlation function. That is equivalently a 
measure of the variability of the signal in the data spec- 
trum. 



3.3. Application to this work: breaking the degeneracy 



When working with low SNR spectra over a reduced 
wavelength range, the cross-correlation functions often 
contain several peaks at different redshifts of similar am- 
phtude. Since the noise reduces the amphtude of the 
peaks of the cross-correlation by a factor directly linked 
to the SNR in the data spectrum, the value of the cross- 
correlation alone is not enough to decide between several 
candidate redshifts and candidate templates. The addi- 
tional information we need to decide on the best match 
is provided by the shape of the cross-correlation func- 
tion away from the peak. A given template is a good 
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Fig. 6. — Spectra of the 25 sources extracted from FIELD 2. In shaded gray we plot the Icr-dcviation for each spectrum. 



match for a given data spectrum if, under the same ob- 
servational constrains (SNR and spectral coverage in this 
case), its "auto-correlation function" provides a good fit 
to the cross-correlation function over the whole range of 
redshifts. 

Each maximum of the cross-correlation function yields 
a candidate solution consisting of a redshift-template 
combination. We use the template and redshift to sim- 
ulate the observed spectrum and compute Pf(^^ j t(-^)' 
that we call an "auto-correlation function" as above: 

Cov{f{z^a^),T) 

Pf{z^^^),TKz) = — — (4) 

We use Fig[TT]to pick the amount of noise to add to the 
template to reach the SNR corresponding to the cross- 
correlation amplitude of the candidate solution. The 
comparison between the cross-correlation and the "auto- 
correlation" tells us whether the choice of z and tem- 
plate recover the shape of the correlation at lags other 



than zero, or equivalently redshifts other than the best- 
fit z. We compute the chi-square between Pf(^zmax),T (^) 
and ps.t{z) for each {z,T) candidate-solution. The min- 
imum chi-square, finally, gives us the best matching re- 
sult which we call the IRS redshift zms- We illustrate 
in Fig [12] the two-parameter nature of the fitting, which 
looks simultaneously for the best-fit redshift and best-fit 
template. This method was automatically applied to all 
extracted SUUSS spectra, without priors. 



4. RESULTS 



4.1. 



Redshift Distribution 

The previous method is successful in determining the 
redshift for 28 sources (62% of the sample). For the 
remaining sources, the lack of redshift determination is 
mainly due to low SNR and/or the absence of any Mid- 
IR feature in our spectral window. Results are presented 
in Table |4] We note that we have a comparable rate of 
successful redshift determination in Field 2 at 14-2 1/im 
(17/25 sources or 68%) and in Field 1 at 20-35/L(m (11/20 
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Fig. 7. — On the left, maximum SNR achieved in each of our spectra versus the observed wavelength at vfhich it peaks. On the right, flux 
density fi, of the SUUSS sources where their spectral SNRs peak. In both plots, the bars in the wavelength direction relates to the contracted 
total spectral coverage over which each spectrum SNR is higher than 2. Empty symbols are sources for which MIPS /i/(24^m) > lOOfiJy, 
and filled symbols are sources for which MIPS fv{2ifiTn) < lOOfiJy as well as sources for which we do not have a MIPS 24fim flux. 

termine mid-IR spectroscopic redshifts as accurate as 1% 
of (1 + z). Where an optical spectroscopic redshift is 
not available, and the MIPS flux density i s higher than 
SSjjJy, we use photometric redshift from |Caputi et al.| 
(20071. The seven such sources in our sample yield a 
< Az/(1 + z) > of 1.5 X 10~^ with a standard deviation 
of 6.6 X 10~^. This dispersion is higher than expected 
from the ph otometric redshift characterization in |Caputi] 




et al.|(2007l. 

Ut the 17 sources without an IRS redshift 5 have an 
optical spectroscopic redshift and 12 have no redshift at 
all. As stated in Sect. 2.5, we do detect all sources with 
MIPS S24 > 80/i Jy in our sample thus reaching compara- 



ble depth as the sample used in Caputi et al. 



Fig. 8.— Ratio of IRS derived and MIPS 2Atim. fluxes 

against MIPS 24/im flux for the SUUSS sources from Field 1 (which 
overlap MIPS 24/im band). The dotted line is the 1-to-l reference. 
The red dashed li ne is the linear fit between 524, /hs E'nd S24,mips 
presented in Sect \2.5\ 

sources or 55%). We flnd a median redshift of 1.05 
slightly larger than the median redshift of 0.935 fo und 
for the 24/ i m flux -hmited sample {S24 > 80// Jy) by Ca- 



(,2007|) 



which 



was selected to have only sources with A'24 > SO/^Jy. Fig- 
ure [13] show^_ou££edshiftdiste and we compare it 
to that of iCaputi et al.l (|2007|) . 



puti et al.| ( |2007[ ). Interestingly, mean redshifts for each 
held are quite different and are equal to 1.25 and 0.92 in 
Field 1 (LLl) and Field 2 (LL2) respectively. We identify 
more higher redshifts sources in Field 1 (4 / 20 sources 
with z > 2) than in Field 2 (1 / 25 sources with z > 2) 
which is probably due to the different rest-frame wave- 
length ranges observed by LLl (6.7 to 11. 3 fim at z = 2) 
and LL2 (4.7 to 6.7/im at z = 2). The most distant 
galaxy is at z ~ 2.2. 

We compare the IRS spectroscopic redshift of ou r sam- 
ple w i th the optical spectroscopic redshifts from Wirth' 
et al.| ( |2004] ) (13 sources, see Table [I] and Fig[l3|. The 
mean of Az/(1 -I- z) is equal to 4.3 x 10"'^ and its stan- 
dard deviation is of 1.0 x 10^^. This shows that, even 
with very noisy IRS spectra, our method is able to de- 



4.2. Spectral Type 

As presented above, the cross-correlation method pro- 
vides spectroscopic redshifts along with an indication of 
the most similar spectral type. We have classified the 
sources in our sample by the best-fit template found 
for each. We divide the sources among four roughly 
defined populations: those distinguished by prominent 
aromatic emission features (14/45 sources), those charac- 
terized only by the presence of silicate absorption (4/45 
sources), an intermediate population gathering mixed 
signatures from weak PAH and/or silicate absorption 
(8/45 sources), and those for which we could not find 
a redshift, and thus a spectral type (17/45 sources) ei- 
ther because of very low SNR or because no features fall 
into the spectral band covered by this survey (14-21/im 
or 20-35/xm). An absence of features in the spectrum 
can result from a high redshift (z > 3.5 in LLl or z > 2 
in LL2) or the intrinsic properties of a continuum dom- 
inated source usually associated with a dominant AGN. 
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(1) IRS-extracted 24 micron 




Fig. 9. — On the left, an IRS-computed 24/im map for Field 1. The contour levels were computed in logarithmic scale in order to 
highlight faint sources. The first contours correspond to 0.3 rr. Right, the contours extracted from the left panel are overplotted on the 
GOODS-North MIPS 24fim observations (Dickinson et al., 2003). 



Two sources for which we have been able to get a redshift 
are not part of any of the subsets described above. We 
will discuss these two in greater detail in the next section. 
It has been previously shown that H2 lines can sometimes 
be detected even in low-resolution IRS spectroscopy (e.g. 
Armus et al.||2006[ iDasyra et al.||2009| . We do not how- 
ever detect any reliable H2 emission m our spectra. We 
show the complete spectral type statistics of the SUUSS 
sample for each field in Table 5j 

A majority of the identifica matches (14/28 sources) 
in our sample present some PAH emission features. The 
spectral ranges available can not cover all the PAH bands 
at the same time, we refer to it as partial PAH emission. 
In spectra with limited mid-IR spectral coverage (in this 
case 14-21/xm or 20-35/im), there is a degeneracy between 
a low SNR 7.7/iTO PAH feature on the one hand and sil- 
icate absorption that creates a false "continuum bump" 
around 8/im in the other. The easiest way to overcome 
this would be to ascertain the existence or absence of an- 
other PAH feature (e.g. 6.2 or 11.3/j,to) not affected by 
the presence of SiO absorption. This is however almost 
impossible here due to the rcdshifts of our sources that 
tends to throw 6.2 or ll.S/im rest frame wavelength out 
of the observed spectral window. 

Figure 14 shows spectra in restframe (using zms only) 
together with the best template. Templates have been 
normalized to the mean fluxes of our sources computed 
over the whole available spectral range. Provided a 
spectrum with high enough SNR (empirically > 3) and 
enough dynamic range (measured as > 1, cf Eq. 2 
and sect. 3.2 1, the effectiveness of the cross-correlation 
method will only be limited by the diversity of mid-IR 
SED shapes covered by our templates. These condi- 
tions lead to very good agreement as seen, for exam- 
ple, for sources SUUSS 20 or SUUSS 3 in FigfTH Sharp 



bidden lines (e.g. [Nelll], [OIV]). On the other hand, 
wider cross-correlation maxima (usually ranging from 
high anti-correlation to high correlation) reveal a more 
diffuse identification based on the global shape of the 
spectrum ( "lower signal frequencies" ) that we can see 
with low SNR 7.7/im emission and SiO absorption (e.g. 
SUUSS 3, 41 or 42). 
Analysing results fr om Table [4| and Table [5) we see that 



cross-correlation peaks at the good redshift denote higher 
SED-feature "frequencies" (e.g. SUUSS 20, 30 or 25) 
such as 11.3/im and 12.7/xm PAHs and/or atomic for- 



PAH templates from Smith et al. ( |2007b[ ) provide most 
of the identifications (14 sources). I'hose templates show 
PAH emission (6.2, 7.7 and/or 8.6//m features, 12. 8/im 
emission line, etc ), a "higher- frequency" signal (com- 
pared to continuum and silicate absorption) that provide 
sharp cross-correlation peaks and thus more accurate and 
reliable redshift identification. A few other templates 
also match PAH emission features such as NGC6240 (2 
sources) or UGC5101 (1 source) with differences being 
mainly on the continuum and the 11.3/7.7 ratio. Other 
identifications mostly rely on broad spectral behavior 
and thus preferentially match templates such as Arp 220, 
Mrk231, Mrk463 or IRAS 15250 (6 sources total). A few 
sources remain with in-between signal properties such as 
Mrk273 (2 sources) . Overall, the cross-correlation analy- 
sis provides strong characterization of mid-infrared emis- 
sion features in our sources, that we propose to use as 
a first order diagnostic to distinguish between starburst 
and AGN as power agent for these sources. For the fea- 
tureless spectra the cross-correlation method does not 
even yield a slope, and therefore generates limited infor- 
mation about the nature of the sources. We computed 
restframe composites for each group of sources with a 
redshift and similar spectral type. The increased SNR 
and spectral coverage are expected to allow us to deter- 
mine additional global properties for those populations 
(cf Sect[5l]). 

No spectral type information could be extracted for 
the 5 sources for which we have an optical spectroscopic 
redshift but no IRS redshift. At least 2 of those have 
a mid-IR spectrum "polluted" by a spatially close and 
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Fig. 10. — Full IRS Low resolution spectra (SL and LL modules) of all the templates used in the Cross-correlation analysis. The 5 PAH 
templates are presented in Smith et al. 2007. IRAS sources (05189, 08572, 12112, 14348, 15250, 22491), Arp220, NGC6240, UGC5101 and 
Mrk 231, 273, 463 and 1014 are part of the IRAS Bright Sample. Their IRS Mid-IR spectra are discussed in Armus et al (2006). The 2 
QSO spectra (3cl20 and PG1612+261) are from P. Ogle (private communication). NGC1569 spectrum is presented in |Wu et al.| | [2006| l. 



brighter source also present in our sample. While we 
have been able to separate spectra of confused sources in 
the previously described redshift identifications, in these 
two cases the "polluting" source was too strong in com- 
parison and contributions could not be spatially or spec- 
trally separated. For one of them (SUUSS 19) there is 
a hint of PAH emission (11.3/im and 12.7/zm) in agree- 
ment with the existing optical spectroscopic redshift of 
the counterpart of z — 1.013. 



4.3. Two special line-sources: fSiv] and [Oiv] lines 



Among the detected sources, two (SUUSS 9 and 17) 
have only one prominent emission line each. In both 
cases, the redshift determination method had to be ap- 
plied more carefully. 

We observe SUUSS 9 at position 189.210 +62.2869 
with a spectrum showing only one emission line at 
32.40/^m (see Fig 15 1 with a total flux of '--^ 8 x 
10~i9 W m~^. A spurious origin for this feature was 
ruled out by several means, one of which was to build 
8 separate cubes (one from each AOR). The feature 
was found in 5 out of the 8 cubes, which is consis- 
tent with the lower SNR when using only one eighth 
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Fig. 11. — Value of a peak in a cross-correlation function as a 
function of the SNR in the data. We used a PAH template from 
Smith et al. (2007) redshifted at ^ = 1. The modified template 
(with noise added) was cross-correlated with the original PAH tem- 
plate and the value of the peak (at z = 1 in the cross-correlation 
function) was tracked as a function of the noise. At each SNR, 30 
realisations of the noise have been computed. The red curves show 
the mean it la of the distribution at each SNR value. 



of the data. The reliability of this source was further 
verified by inspection of its spatial profile. We inte- 
grated fluxes in a Hanning-window centered on 32.4/LtTO 
and extracted the source PSF which matches that of 
brighter sources in our sample. We address the problem 
of assigning a redshift based on a single line by relying 
on our knowledge of typical mid-IR spectra of known 
sources, and by using surveys at other wavelengths in 
the same part of the sky f or additional evidence. There 
are three GOODS sources ( Chary|2007 l sufficiently close 
to be considered as potential counterparts, as reported 
in NED (NASA/IPAC Extragalactic Database, Decem- 
ber 2008): GOODS J123650. 87+621712. 8 is 3.4" away 
from our estimated position of SUSS9, and at a spectro- 
scopic redshift of 2.133; GOODS J123650.22-I-621718.4 
is 5.5" away from SUSS9, and at z = 0.51283; GOODS 
J123649.44-h621712.3 is 6.7" away from SUSS9, and at a 
p hotometric red shift of 0.04. 

( |2007b| and |Dale et al.| (|2006| have 



Smith et al. 



reported the detail ed spectroscopi c cont ent ot SINGS 
galaxies; similarly, Armus et al. (20071 reported on 
nearby ULIRGs. We will examine which of the most 
prominent lines seen in those systems might be associ- 
ated with this one line, in view of the mid-IR spectral 
properties of those samples, and of the visible sources 
positionally associated with SUUSS9. We should clarify 
however that mid-IR spectra are still yielding many sur- 
prises, as in the case of radio galaxies (Ogle et al 2009, 
in prep), so that plausibility arguments presented here 
cannot yet be quantified as to their likelihood. 

Molecular hydrogen is known to dominate the mid- 
IR spectra of certain object s, as in th e case of Stephan's 
Quintet intergalactic shock ( Appleton et al. 2006), or the 
Radio Galaxy system 3C 3261 Ogle eral.i2007| ) . One pos- 
sibility is that this could be the ti2 S(l) line at z — 0.902; 
one would then expect the S(2) line at 23.35/im and the 
[Neil] line at 24.36/im, neither one of which is detected. 
In Stephan's Quintet, S(l) is the most luminous line, 
and each of [Neii] and S(2) carries one third or more of 
the luminosity in S(l), so their absence in SUUSS 9 is 
a significant argument against this possibility, as is the 
absence of any nearby visible sources at redshifts close 



to 0.9. Similarly, this could be the H2 S(0) 28/zm line 
at z = 0.148, with the S(l) line expected at 19.55/im, 
outside the survey spectral coverage. The lack of visible 
sources in the redshift vicinity of 0.148 argues against 
this possibility. Finally, if we assume this is the H2 S(2) 
line at z = 1.639, then the S(3) would be expected at 
25.5/j,m, and it is not detected. The weak signal peak- 
ing around 16.97/im (restframe) is significantly removed 
from the expected wavelength. This non-detection is a 
strong argument given the typical weakness of the S(2) 
line, so we conclude this is an unlikely assignment. 

The most prominent fine-structure lines that might be 
associated with the detection at 32. 40 /irn are: [Siv] at 
z = 2.083, [Neil] at z = 1.529, [Nev] at z = 1.262, 
[Neiii] at z = 1.083, [Siii] at z = 0.731, [Oiv] at 
z = 0.251, or [Fell] at z = 0.247. Among all these pos- 
sibilities, we favor the [Siv] assignment, for two main 
reasons, one being the closest of the visible sources has 
a very similar redshift at z = 2.133, and the other being 
the existence of a similar mid -IR spectrum fr om the local 
Universe, namely NGC1569 ( |Wu et al.|[2006| . The other 
potential assignments are less likely because the observed 
spectrum does not contain expected lines or Aromatic 
features, and none correspond to any of the redshifts of 
nearby sources. For instance, [Neii] is quite unlikely to 
appear without Aromatic features, and we would have 
detected the 11.3/im band at 28.6/zm for z = 1.529. How- 
ever, the radio galaxy 3C317 does display a high ratio of 
[Neil] to Aromatic Features (Ogle et al 2009, in prep). 
Similarly, [Ne v] is generally much weaker than [Ne 11] , 
and the latter would have been detected at 26.7/um for 
z = 1.083. [Neiii] on the other hand can be much more 
luminous than [Ne li] , thus appearing as the only line in 
a spectrum, but only in a small minority of cases. Dale 
et al (2009) find that the [Ne ill] flux exceeds twice the 
[Ne 11] flux in less than 10% of the cases for a variety of 
SINGS objects. While it is possible that the line we de- 
tect is [Ne III] , and that [Ne 11] has escaped detection at 
28.99/^m for z = 1.262, this redshift does not correspond 
to any objects in the fleld. The same argument applies to 
[Siii], since [Neii] would have been detected at 22.19/im 
for z = 0.731, and it is rare for [S ill] to exceed twice the 
flux of [Nell]. The lines of [Oiv] or [Fell] near z = 0.25 
are again unlikely candidates, since the [Siii] line is ex- 
pected to be at least comparable in flux and detectable 
at 23.4/im. The exception would be that [O iv] might 
exceed [Siii] in AGN-dominated sources. However, the 
lack of a detectable continuum and the lack of a similar 
redshift in this vicinity are both arguments against this 
identiflcation. 

The close association of SUUSS9 with GOODS 
J123650.87-I-621712.8 both spatially and in redshift is 
suggestive of these being the same object, but not con- 
clusive. It is just as likely that they are two neighboring 
sources rather than the same source. The redshift dif- 
ference amounts to ~ 15,000 km/s, and corre sponds to 
1.7cr assuming the uncertainty on the Reddy et al. (2006| 
redshift is sufficiently lower than ours. Moreover, inter- 
preting SUUSS9 as a NGC 1569 analog implies that it is a 
Wolf-Rayet galaxy, whereas the GOODS source does not 
displa y the Wolf-Rayet charact eristics in its optical spec- 
trum ( Brinchmann et ar]|2008[ ). Additional data will be 
required to clarify the relation between these two sources, 
which is beyond the scope of this paper. While we can- 
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Fig. 12. — We represent the cross-correlation between 3 different spectra and 5 different templates witfi very diverse mid-lR spectral 
signatures. In each case we show the cross-correlation function (black) and the auto-correlation function (red) for the best matching redshift 
(red square) with that specific template. The best candidate solution (redshift and template) for each source is shown as a gray shaded 
plot. 
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Fig. 13. — Left, IRS spectroscopic redshifts (zjjig) from our study against existing spectroscopic (TKRS, Wirth et al., 2004) and 
photometric (Caputi et al., 2007) redshifts (z^^/)- Right, redshift distributions. The dashed-gray histogram show the total redshift 
distribution for our sources {zjug supplemented with Zj-e/ when z^^s couldn't be obtained). The shaded histogram shows the sources for 
which we had only z^-^f. The dashed open histogram is Caputi et al. (2006) sample distribution. 
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Fig. 14. — On the left, our sources spectra (black) with best match template (red) overplotted. The template is scaled to the mean flux 
of the spectrum. On the right, the related cross-correlation function (in black) and auto-correlation function (in red). Results for the first 
4 sources are shown here as an example, plots for the whole sample can be found at the end of this document. 
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Fig. 15.— Spectrum of the [SIV] "line-source" (SUUSS 9) in 
restframe. Icr-deviation is shown as the gray shading. The cross- 
correlation gave us a redshift of z = 2.08 for this source. In dashed 
red we show the spectrum of NGC1569 for which we obtained the 
best cross-correlation match. It is scaled to our source 10.5/^m 
restframe flux density. 



not attach a formal confidence level to this statement, 
we do propose that the interpretation of the SUUSS 9 
line as [Siv] at z ~ 2.1 is the most plausible postulate, 
followed by the less likely interpretation as the H2 S(0) 
Une at z = 0.15. If SUUSS 9 is indeed detected in [Siv], 
then its [Siv] line luminosity would be 8.3 x 10^^ W, 
whereas the dwarf galaxy NGC1569 emits 1.5 x lO'^^ W. 
While this is a dramatic difference, the [Siv] luminosity 
estimated for SUUSS9 is not extraordinary, in the sense 
that the quasar PG1612-f 261 has a comparable luminos- 
ity of 1.4 X 10^^ W. The interpretation of the SUUSS 
9 line as H2 S(0) a,t z = 0.15 is less demanding ener- 
getically, since it requires scaling up the emission of the 
intergalactic shock in Stephan's Quintet by a factor 11 
only, to a line luminosity of 1.5 .lO"^'^ W. 

The second line-sources (SUUSS 17) possesses a rel- 
atively well determined mid-IR continuu m o f about 
SOO^iJy (> 
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This 



6a) as can be seen in Fig 
source shows one very well defined emission line at 
A ~ 33.7/iTO (observed wavelength) with a total flux of 
~ 6 X 10^^^ W m^^, and no other significant spectral 
feature over the LLl wavelength range (20 to 35/xm). 
The two absorption-like features at 26 and 28.5/im in 
the spectrum are unlikely to be real and are more likely 
originating from a combined effect of noise and confu- 
sion as this source sit in a denser region of our survey 
with 3 surrounding sources as close as 15". The line is 
well resolved by the IRS, with an estimated intrinsic full 
width at half max equivalent to ~ 4000 ± 1400 km/s. 
NED, consulted in December 2008, reports three objects 
sufficiently close to be considered potential counterparts: 
GOODS J123658.45-I-621637.3 at z = 0.2993, and 1.5" 
away, GOODS J123658.82-H621638.1 at z = 0.29863, 
and 3.2" away, and GOODS J123658.09-H6216 39.4 at 
z = 1 . 01734 and 6.2" away; all redshifts are from Wirth| 
et al. ( 2004 ) . The second of these is assigned a Spitzer 



24/im detection at 264/iJ?/, making it the most credible 
counterpart to SUUSS 17. 

The cross-correlation analysis on this source yields sev- 
eral candidate solutions for different templates and does 
not allow us to secure a reliable redshift. However, using 
the z = 0.3 redshift as prior we find the best match- 
ing template to be the spectrum of pgl612-|-261 (QSO, 



Shi et al.|[2007| with an IRS determined redshift of 0.3. 
This redshift implies that the strong emission line in the 
IRS (LLl) spectrum is either [Oiv] (25.89/xm) or [Fell] 
(25.99fnn). Moreover, at z = 0.3, the IRAC 8^m band 
falls on top of the 6.2fim feature and would be expected 
to cause an excess compared to the shorter wavelength 
IRAC bands. Such an excess is evident in the inset in Fig 



16 and we tentatively interpret it as the presence of PAH 
emission. The difficulty with this interpretation however 
is that [Oiv] is ob served with < 1000 km/s width in 



known QSOs (e.g. Dasyra et al.| |2008|, and the line 



width observed for SUUKIS 17, typical of broad-line re- 
gion emission, is not expected for high-ionization species 
like [O iv] . It is also difficult to explain the large observed 
line width as resulting from the combination of [Fe 11] and 
[O iv] emission, for the separation between the two lines 
is only O.l/im, or 900 km/s. One might invoke the pos- 
sibility of a broadened [Fe 11] rather than the [O iv] , but 
one would then expect the 18.7 fj,m line of [Siil] to be 
at least comparable, and typically several times stronger 
than [Fell] (Dale et al 2009, in press). There is no ev- 
idence of the [S III] line in Fig 16 which argues against 
this possibility. The 17.94/iTO Ime of [Fell] is typically 
weaker than the 25.99/im line, so we would not expect 
it to appear. Another strong argument against SUUSS 
17 displaying a broadened [Fell] line is that this this line 
is ra rely seen as th e dominant emission line from QSOs 
(e.g. Dasyra et al.| [2008), and the lack of a strong [Siii] 
line i s more consistent with the observed line being [O iv] 



(e.g. lOgle et al.l|2006|) 
■ ^ ■bii:' 



Another possibility is that the 33.7/im emission is an 
Aromatic feature rather than a broadened line. We find 
this improbable, since other PAH features would be ex- 
pected in the same spectral range: for example if we 
assume this emission feature is the 11.3/im PAH (at 
z = 1.98) then emission at 7.7/im restframe should be 
observed at A = 23/xm. While the overwhelming ma- 
jority of sources follow this expectation, Leipski et al 
(2009, submitted) show a clear counter-example in M84 
(3C272.1), a FR-I radio galaxy with an IRS spectrum 
displaying the 11.3/im feature but none of the features 
at 6 to 9/im. However, given the strong evidence of an 
independent z — 0.3 determination for a spatially coin- 
cident source with the matching continuum flux density, 
we conclude that this possibility of z = 1.98 is much less 
likely than the z — 0.3 identification. 

With the [O iv] aX z — 0.3 identification, the large line 
width suggests the presence of a strong QSO. GOODS 
J123658.82-H621638.1 i s also detected with Chandra 
(Alexander et al. 2003). This is a soft X-ray e mitting 



source and falls in the "normal galaxy" category of |Bauer 



et al.| (2004) cl assification . In a more recent analysis by 



jcorgakakis et al. (2007), the ratio of x-ray to infrared 
luminosity also places this source in the middle of the 
range for star-forming galaxies. Note also that promi- 



nent Aromatics at 6.2fj,m (Fig 16 inset) also argues for 
a strong star formation contribution to the luminosity of 
this source. This apparent contradiction with the mea- 
sured line width in the spectrum may be the result of high 
dust attenuation affecting the x-ray luminosity, or sim- 
ply due to multiple sources contributing to the IRS spec- 
trum. The latter would not be surprising as the GOODS 
catalog contains two sources aX z ^ 0.3 within 3.2" of 
SUUSS 17. Here again, the additional work needed to 
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Fig. 16.— Spectrum of the [OIV] "line-source" (SUUSS 17) in 
restframe. Itr-deviation is shown as the gray shading. The green 
square is the MIPS 24/im counterpart flux in restframe. The cross- 
correlation, in agreement with ancillary data, gave us a redshift of 
z = 0.3 for this source. IRAC data points (inset) tend to agree with 
the presence of PAH emission at 6.2 fim. In dashed red we plotted 
the spectrum of the QSO PG1612-I-261 for which we obtained the 
best cross-correlation match. It is scaled to our source 25.9^m 
restframe flux density. 

resolve these ambiguities is beyond the scope of this pa- 
per. 

5. DISCUSSION 

This survey offers us a way to characterize the Mid- 
IR population of galaxies in a small region of the sky 
down to very deep detection levels. While any survey 
samples the Universe imperfectly, this survey avoids the 
bias of a narrow bandpass sampling complex SEDs. No 
bias is expected from the detection method used as it is 
only related to the integrated fluxes of the sources in the 
field and does not rely on any specific spectral feature 
detection. 

We investigate the possibility that the analysis method 
presented here favors sources with specific types and red- 
shifts in our field, potentially resulting in an a posteri- 
ori bias in the results. As explained in Sect. 3.1, the 
cross-correlation method is sensitive to the variance of 
the signal in the spectral window with respect to the 
variance of the noise. This means that the method gives 
better results when applied to spectra that show signifi- 
cant structure such as strong emission lines. Aromatic 
Features or silicate absorption. Considering that the 
spectral window in our sample is limited to either 14- 
20/im or 20-35/im this could lead to inhomogeneities or 
biases over redshift and type of sources for which the 
cross-correlation method was successful. If this effect 
is present, the only efficient way to remove it is to add 
multi-wavelength data (IRAC fluxes, near-IR and optical 
bands ... ) to the spectra used in the cross-correlation 
algorithm thus lowering the probability of missing signif- 
icant spectral information. 

We discuss below the luminosities of the sources de- 
tected in this survey, then their mid-IR spectral charac- 
ter, and place them in the context of other surveys. We 
then contrast this unbiased spectral survey with single- 
band continuum surveys. 

5.1. IR luminosities and spectral types 

Since we only have mid-IR data, we can not compute 
accurate total IR-luminosities. We could use the best 
template for each source but some of the templates have 



no far-IR data and 42% of the sources are not identi- 
fied. We estimate IR bolometric luminosities from MIPS 
24/im flux measurements using a code made available 
online by R. Char y. The code sele c ts th e most ap pro- 
priated SED from Chary fc Elbai| (f2(501 ) and [Dakfc] 

XJ^l^,, /onno\ rTTTTT^n TTT^ iU J 0/1,,™ ^ 



Helou (2002 ) models based on the 24/im flux, fits the 
SEDs and computes the total 8-lOOOum IR luminosities, 
then averages them together. Figure 17 shows the source 
luminosities as a function of redshift. The luminosities 
range from 3.2 x IO^Lq to 4.1 x lO^^L© with a mean of 
5.2 X 1O^^L0. Even the 4 sources at z > 2 have an im- 
pressively low mean luminosity of 7.5 x IO^^Lq. When 
compared to existing IRS spectroscopic studies of high- 
redshift galaxies (Houck et al. 2005, Pope et al. 2006, 
Yan et al. 2007), our sample proves to be one of the 
deepest in this range of work (see Fig 
We selected all sources at 0.8 < z < 
at 1.7 < z < 2.2 (4 sp ectra, the 
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1 spectra) and 
"line-source" SUUSS 9 
discussed in Sect |4.3| was discarded from this sample) 
and computed corresponding composites (Fig 18 1. Each 
spectrum was scaled to its sample mean luminosity. The 
composite at z ^ 1 is dominated by Aromatic Features. 
We detect the 7.7, 8.6, 11.3 emission features as well 
as the [Neil] 12.8 fim line. We see generally agreement 
with the AGN-free starburst template from Brandl et al. 
(2006). Due to lack of spectral coverage in this z ^ 1 
sample, we have no data at A greater than 17/L(m. 

In the composite at z '--^ 2 we can clearly see 6.2, 7.7 
and 8.6 /im PAHs. In comparison to z ~ 1, this sample 
of sources is half an order of magnitude more luminous. 
We also can not confirm or rule out the presence of any 
emission at 11.3 /im, although this feature is very often 
detected across the entire sample of this study. This may 
be due to the small number of sources entering the z ~ 2 
composite together with increased noise at this restframe 
wavelength. The shape and relative strengths of the 7.7 
and 8.6 fiia features is however conserved between z ~ 1 
and z 2. The 4 sources that enter this composite all 
show a clear stellar bump in the IRAC channels which is 
expected for starburst galaxies betwee n z = 1.7 (bump 2 



and z — 2.2 bump 3). As seen by Farrah et al. (2008 1 



Tared 



in their sample of 32 high-redshift ultralummous m: 
galaxies selected as bump 2 sources (mainly), composite 
spectra of such sources are best described by a st arburst 
component, namely the template produced by [Brandl 



et al.|(|200 6l 
rah et 
either. 



al. 



Little silicate absorption is seen in the Far 
sample and none is seen in our composites 
The same observation is made on the population 
of submilli meter-selected high redshift i nfrare d ga laxies 



presented in Menendez-Delmestre et al. ( 2007 ) and Pope 
et al. ( 2008 1 where a starburst scenario ('Brandl et al 



2005 template or M82) with low-AGN contributio n"^ 
ter fits the co niposite spectra. P revious work by Houck| 
et al. (2005 1 or |Sajina et al. ( |2007 1 find prominent and/ or 
non-negfigible SlU absorption m their sources (55% and 
33% of the samples respectively) . Our z 2 composite 
is surprisingly quite similar to the com posite of z 2 



sources presented in Sajina et al. ( 2007 1 with the excep- 
tion of the slight SiU absorption. Overall, we observe 
mid-IR composite spectra generally similar to those de- 
rived from other samples of high-redshift (1 < z < 2.5) 
infrared galaxies with very different selection criteria and 
bolometric luminosities, and generally favoring a star- 
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the Chand ra Deep Field North 2Ms catalog (Alexander 



Fig. 17. — IR bolometric luminosities computed for our sample 
from MI PS 24/jm fluxes using a cod e provided by R. Char y imple- 
menting ICharTX^jbazl |[20^ and [Dale fc Helou' f 2002|) models. 
Uncertainties are ot 5U7c. We oyerplot IK luminosities of several 
other IRS surveys of high-redshifts galaxies. Filled dots are SU- 
USS sources from Field 1 and empty dots are SUUSS sources from 
Field 2, yellow filled triangles are 17 ULIRGS from the IRS GTO 

eouck ct al. 20051, blue open triangles are 47 ULIRGs from the 
'LS ( Yan ct al. 2U07), green diamonds are bump 2/3 sources from 
[Farrah ct al. (2U0S;), orange ope n squares are 13 SMGs from the 
ilUi'-Nortli SGUBA Super-map | Pop e et al.l2008|, purple crosses 
are 5 SM Gs from Mcncndcz-Dclmcstrc ct al. (20071), and red stars 
are from [Teplitz ct al. (2007 ). I'hc dotted line is the luminosity 
limit of the bUUbb survey (jl{24:fim) = A7fj.Jy). 

burst character. 

5.2. Characterizing the population (AGN vs starhurst ?) 

Table [5] summarizes the statistics of spectral types 
assigned for objects detected in this survey. We as- 
signed each spectral template used to one of three classes, 
namely strong aromatics, silicate absorption and mixed 
signatures. The counts of objects falling in each category 
are shown in columns 3, 5 and 4 respectively; column 7 
accounts for sources which could not be matched using 
the spectral corre lation method. The two sources dis- 
cussed in Sect 4.3 were gathered in a fourth class (column 
6). Among the sources with accepted matches, about half 
(14 out of 28) are dominated by the aromatic features 
and another third (8 out of 28) have a mixed character. 
About 15% (4 out of 26) have substantial silicate absorp- 
tion, indicative of high optical depth, against a dominant 
continuum suggesting a major AGN contribution. Out of 
the total sample of 45 sources, 17 could not be matched to 
any templates in our library. This might be due to their 
being dominated by a featureless continuum, or to their 
low signal-to-noise ratio. If the former is the dominant 
cause, then these 17 might be primarily AGN-dominated 
sources, and the statistics among matched spectra are bi- 
ased towards starburst signatures. 

In order to rule out a bias in favor of starburst galax- 
ies among matched spectra, we used additional diagnos- 
tics based on IRAC and MIPS color-color plots (Fig[l9l. 
They can be used to separate AGN-dom inated sources 



up to z ~ 2 (e.g. Donley et al. 2008). According to 



these diagnostics, no clear AGlN-dommated system is de 
tected in our total sample of 45 sources, which is consis- 
tent with the redshift and spectral-type determination 
among sources matching templates. The sources in our 
sample therefore seem to be mostly dominated by star- 
formation, and this result is not a bias due to the spectral 
matching technique. 
We have also used X-ray data for this region from 



et al. 2003 1 to look for X-ray counterparts to our sources 
We found four matches inside a 3 arcsecond matching 
radius . Fol lowing the classification proposed by Baue^ 
et al. (12004), two of them are potential X-ray obscurecT 
AGNs (SU USS 31 and 41), and one is a starburst galaxy 
clearly identified as such with our dataset (SUUSS 20). 
The last X-ray counterpart is the [O iv] line-source pre- 
viously described in Sect. 3.4. Three additionnal coun- 
terparts were identified using the supplementary cata- 
log; these are sources that were not detected with a high 
enough significance to be included in the main 2Ms X- 
ray catalog. They are fainter sources but have an optical 
counterpart which makes them highly likely to be real X- 
ray sources and potential candidates for faint AGN. One 
of them presents PAH emission with possible silicate ab- 
sorption, one is clearly continuum dominated (with 17/im 
PAH emission detected) , and the last one has deeper sil- 
icate absorption which give us two other possible AGN 
identifications. Tentatively, we estimate an AGN frac- 
tion of 4 out of 45 (9%) of our sample, based on X-ray 
detections. Given the sample size, this is consistent with 
the spectral typing results sum marized above . 
The IRS GTO observations of|Houck et al.|(2005) un- 



veiled a population of sources dominated by strongly ob- 
scured (SiO absorption) AGNs, so me presen ting PAH 



signatures even though faint. The |Yan et aL (|2007j) 



sample of 52 ULIRGs has more diverse properties with 
roughly a third of the sources presenting strong PAHs, 
another third showing strong silicate absorption and the 
rest being continuum sources sometimes with weak fea- 
tures. These sources show large AGN contributions (80% 
of the sample has some level of AGN continuum signa- 
ture) with a significant fract ion of sources having a star 



burst signature as well (Yan et al. 2007 Sajina et al 



2007| ) . Those two samples have very ditterent AGN trac- 



tions and SiO absorption depths in comparison to our 
sample. However, those samples also derive from quite 
different selection criteria and flux densities. We will ad- 
dress below the relation between the AGN fraction dif- 
ferences and the sample selection. 
Previous studies have shown submillimeter galaxies at 



('Egami 
2006). 



et al. 





2004 




Ivison et al. 


2004 



Pope et al 



|et ar.| (2008) sample may be useful. As we already dis 
cussed, both samples have very similar composite spec- 
tra at 1 < z < 2.5. While X-ray classification tends to 
show the presence of AGN in about 40% of those sources, 
star formation is found to account for more than half the 
bolometric luminosity. An interesting way of com puting 
mid-IR photometric redshifts was proposed in Pope et al. 
(2006) using all four IRAC bands and MIPS 24/^m fKSes 
to give an estimate of the redshift. This method assumes 
a fairly constant SED shape across the sample and thus 
is spectral-type dependent. We produced a similar for- 
mula for our sample by fitting to the redshift obtained 
with the cross-correlation method: 

^^IR = a + & • log(S3.6) + c • log(S5.8) 

+ d ■ log(S8.o) + e • log(S24) (5) 

where the fitting parameters values are 1.7, —2.6, 3.3, 
— 1.0 and 0.1 for a, &, c, d and e respectively. We 
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Fig. 18. — Left, composite spectrum at 2 = 1 it 0.2. Right, composite spectrum at 2: ~ 2. In pray shad ing we plot the la dispersion for 
each composite. In dashed red, we overplot the "starburst average spectrum" from Brandl et al. (2006') scaled to our composites 8.6^m 
and 7.7/i'm luminosities at 2 = 1 and 2 ~ 2 res pectively. We also plot in green average spectrum from Farrah et al. (2008) and in orange 
2 ~ 2 average spectrum from [Pope et al^ l |2008| . Both composites are dominated by PAH emission at 7.7/ir?i and agree very well with the 
starburst template. The ll.Sfim fAU detection in the z ^ 2 composite is not clear. This part of the composite (A > lO/im) is particularly 
affected by the noise. 



find the deviation of the redshift statistic to be of 
(7 {Az/ (1 + z)) = 0.22 which is rather large compa red to 
the di s persio n of a {Az/ (1 + z)) = 0.07 found in Pope 
et al. ( 2006| ) , perhaps explained by the fact that our 
sEUs that enter the fit have a wide variety of spectral 
types. However, no selection could be found that, when 
applied to our sample, reduced the dispersion. It is there- 
fore very likely that there is a larger diversity of SE Ds in 
our sa.mple compared to the 13 SMGs studied by [Pope 
et al. (2008), even though both have very similar average 
mid-lH spectra. 

5.3. Unbiased Spectral Surveys: A Different Probe 

While our survey covers a small solid angle and the 
data are collected over two separate spectral windows, 
it still yields a somewhat different perspective on the 
mid-IR universe at z < 2 and illustrates some of the 
distinguishing characteristics of wide-band spectral sur- 
veys. The basic statistics, summarized in Table l5] point 
to a relative paucity of AGN signature among sUUSS 
sources compared to ot her deep spect r al sur veys such as 
Weedman et al.| ( |2006[ ) or |Yan et al | ( |2007[). The same 
lack of AGIN signatures is observed in IKAC color-color 
diagnostic plots of SUUSS sources, so it is very unlikely 
to reflect a bias introduced by the spectral correlation 
technique. This difference is partly due to additional se- 
lection criteria imposed on the other surveys, partly due 
to the fainter flux densities and luminosities of the SU- 
USS sample, but probably also in part due to the reduced 
selectivity of an unbiased spectral survey. 

The Yan et al and Weedman et al samples were both 
selected for high 24/im/i? ratios, with an additional selec- 
tion for high 24/8/im ratios by Yan et al. These choices 
were aimed at biasing towards high redshifts and lumi- 
nosities, and indeed they achieve that as can be seen 
in Fig [l7j Aromatic Features dominate the spectra of 
no more than a third of those samples. The selection 
for high 24/im/8/im or 24/im/i? favors high obscuration, 
which is mildly anti -correlated with Aromatic signatures 
(Sajina et al. 20071, and might therefore reduce the in- 
cidence of Aromatic Features in those samples. 

On the other hand, the SUUSS sample has a median 
/i/(24/LtTO) ~ 100/iJj/, whereas the other samples are 



brighter by one order of magnitude for Yan et al, and al- 
most 3 orders of magnitude for Weedman et al. The lower 
fluxes lead to lower luminosities, and would generally fa- 
vor star formation powered systems over AGN, translat- 
ing into a preponderance of Aromatic-Feature rich spec- 
tra (Papovich et al 2007). Given the evidence at hand, 
we favor this interpretation as the dominant factor in 
determining the spectral characteristics of this sample. 

The IRAC colors of the SUUSS sample also point to 
a lack of objects residing in the AGN-dominated part 
of that diagnostic plot (Fig 19). The same diagnostic 
is used over all the GOODS- North sources with a MIPS 
24/im flux in between 45 and 100/iJy (Figl20|. Out of 



the 823 sources that fall into this category only 93 (11%) 
are found in the AGN-dominated part of the diagram. 
The effect observed in the SUUSS sample is thus most 
probably due to its lower mean luminosity, which results 
from the sensitivity achieved but also from the choice of 
a particularly dark survey fleld. 

Other factors may affect the character of SUUSS 
sources. One possibility is that spectra with more fea- 
tures, i.e larger 175, are easier to pick out in a spectrally 
dispersed survey than in a traditional single-band sur- 
vey, since a feature would have to coincide with the sur- 
vey fllter to generate a signal. Thus single-band sur- 
veys would favor continuum-dominated sources at all red- 
shifts, bias against those whose absorption features fall 
into the band, and bias for sources whose emission fea- 
tures fall into the band. This may explain the prepon- 
derance of AGN in 24/im selected samples (e.g. Daddi et 
al 2007) . Fig [tI right-hand panel illustrates the lack of 
selectivity in tnis survey, as it shows most sources peak 
away from 24/im, and indeed have higher fluxes where 
detected than at 24/im. 

The 24/im survey obviously missed the single-line 
source SUUSS 9, since the line fell outside the 24/im 
band. Since we found one such source in the whole sur- 
vey, the incidence of such sources in general must be of 
the order of 2 — 4%, depending on what we assume for 
the efficiency of our search. Since this source does not 
have a 24/im counterpart, we can only assign an upper 
limit to its luminosity as deflned in Fig |17[ That upper 
limit is roughly 5 x 10^^ Lq, based on the 24/tm detection 
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t lG. Ly. — spitzer c olor-color diagram characterization ot the IKb Ultradeep survey s ample or galaxies. Lett, ACjlN diagnostic proposed 
by [Pope et al.| | |2008| l determined from a sample of 13 SMGs from zr^l to 2.5 (see also |Ivison et ar]|2004[ l. Right, AGN diagnostic using 
IRau colors, explained in Stern et al. (2005). Tlie Broad-Line AGN selection criteria are represented by the dotted line. The lower and 
left-hand axes are scaled to Vega magnitudes whereas the top and right-hand axes refer to AB magnitudes. 

luminosities from about 7 x IO^'^Lq at z = 0.5 to about 
at z ^ 2.0 poorly covered by IRS spectroscopy so 
far. At least 47% of our sample show starburst activity 
as PAH signatures (up to 21/45 sources) and 31% (14/45 
sources) prove to be starburst dominated. The rest of the 
redshift identifications rest on silicate absorption feature 
(11% or 5/45 sources). We tentatively identify only 9% 
of our sample (4/45 sources) as AGN candidates. This 
small fraction of AGN in comparison to previous IRS 
spectroscopic surveys of high-redshifts galaxies is likely 
to originate from the fainter luminosities reached here. 
We find two unusual sources with only one prominent 
emission line detected over the spectral range covered. 
One of them present a strong [OIV] line (25.9/im rest- 
frame) at z = 0.3 on top of a significant continuum 
while the other shows the [SIV] line (10.5/im restframe) 
at z = 2.08 on a very weak continuum dominated by 
instrumental noise. 



[5.8] - [8.0] (Vega) 

Fig. 20.— IRAC color-color diagram for all GOODS-North 
sources with identified MIPS 24/im flux such as 45 < /^(24/im) < 
lOOfiJy. 823 sources fall into this category and only 93 of them 
enter the AGN region of the diagram i.e 11% of the sample. This 
low ratio indicates the importance of the detection limit on the 
ratio of AGN-dominated system over starburst dominated system 
found in flux limited surveys. 

limit. 

The wavelength coverage in this survey was either 14 
to 21/xm or 20 to 35/xm. This limited spectral coverage 
has contributed to the difficulty of identifying a spectral 
match for more than a third of the detections, and to the 
difficulty of identifying the single line detected in SUUSS 
9. Future IR spectral surveys would be more productive 
if they covered at least a full octave in wavelength. 

6. CONCLUSION 

We have obtained IRS spectra for 45 sources of very 
faint IR galaxies and determined redshifts for ~ 60% 
(28/45 sources) of them using a specifically designed 
cross-correlation method. We covered a domain of IR 
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A Spitzer Unbiased Ultradeep Spectroscopic Survey (SUUSS) 



TABLE 1 
Source extraction information. 



ID 




Field 


RAiBs 


DECiRs 


SNR ^ 


Peak SNR ^ 


Peak Wavelength ^ 


AA2a ^ 








degrees 


degrees 




fim 










LLl 


loy . 1 / u 


RO OQQ/I 

uz .zoyi 


lO.O 


^ 7 


32 22 


A R7 




9 


T T 1 


ioy . loo 


(in OQAQ 

DZ.Z0'4:0 


■4:0. y 


10 9 
lU.Z 


Ol .D / 


1 1 A'i 
11. lo 


o u u oo 


Q 

o 


LLl 




RO 9Q7Q 
UZ.ZO I o 


o.u 


7 ^ 


9"^ on 
zo. uu 


9 QO 


o u u oo 




LLl 


1 QQ 1 Q1 


RO 9QQC; 

uz.zooo 


7 7 


A 1 

4:. / 


ol .u / 


^A 


o u u oo 


o 


LLl 


1 en 1 QQ 
loy . lyo 


RO 9QA4 

UZ . ZO'4:4: 


0. 1 


O.O 


9fl 

ZO.oO 


1 77 
1. / 1 


o u vj oo 


u 


LLl 


loy .zuo 


RO 9QA7 
uz .zy^i i 


2 2 


1 ^ 

l.O 


Ol .U ( 


U.IO 


o u VJ oo 


7 


LLl 


loy .zuo 


RO 97fi/i 

UZ . Z i U4: 


7 Q 
/ .y 


^ 9 
o.z 


29 11 


o.uo 


o u u oo 


c 
o 


T T 1 


loy .ziu 


(^9 9S97 
DZ.ZoZ 1 


o.o 


9 
D.Z 


ol .oO 


O.OO 


SUUSS 


gd 


LLl 


189.210 


62.2869 


1.2 


3.2 


33.50 


0.16 


o u u oo 


1 n 

-LU 


LLl 


1 ^iQ 917 


R9 QQ'^R 


o.o 


o.y 


31 12 


1 77 
1. / 1 


SUUSS 


11 


LLl 


189.225 


62.2907 


11.5 


8.7 


22.65 


5.64 


OTTTTCO 




T T 1 


1 on TOO 


Dz.zozy 


27.5 


8.6 


27.84 


7.89 


QTTTTQQ 

oU Uoo 




T T 1 


1 on OQO 

iisy. zoz 


Dz.zozz 


Q Q 
O.O 


7 1 


z4.4o 


/tin 

4.±y 


oU Uoo 


14 


T T 1 


loy.zo I 


oz.z i oo 


19.3 


6.3 


24.07 


7.25 


CT TT TQC 

oU Uoo 


15 


T T 1 


ion O /I Q 
ioy.z4o 


DZ.2DOO 


22.6 


6.2 


33.50 


7.7o 


CT TT TCC 

o U U oo 


Id 


T T 1 

LLL 


1 en o/i /I 
ioy .z44 


Dz.zDo4 


1 n 
iy.o 


O.z 


oo. OU 


( .Uo 


CT TT TCC 

oU Uoo 


1 / 


T T 1 

LLL 


ion o A A 
ioy.z44 


(^O 0'7'71 
DZ.Z ( ( L 


67.6 


15.4 


33.32 


12.56 


CT TT TCC 

o U U oo 


1 Q 

io 


T T 1 

LLL 


1 cn o/i R 
ioy .z4d 


Oz.zoiU 


u.o 


/ .4 


Qn K'T 


0.04 


CT TTTQQ 
o U U oo 


1 o 


T T 1 

LLL 


1 fiQ Qf^l 

ioy .zoi 


Dz.ZDo ( 


Zi .4 


O.l 


Q1 Qc; 
oi.ou 




OTTT TCC 

oU Uoo 




T T 1 

LLL 


1 on oc^ 1 

ioy.zoi 


DZ.Z ( i4 


QQ 


on /I 
zy.4 


z4.4o 


12.00 


OTTT TO C 

oU Uoo 


21 


LL2 


189.259 


OZ.ZOOZ 


14.3 


6.6 


15.28 


3.08 


OTTT TO O 

oU Uoo 


22 


LL2 


"ion o o 
l09.ZDO 


/"'o o cr 1 n 

Dz.zoly 


10.3 


4.3 


20.32 


1.08 


OTTT TOO 

oU Uoo 


23 


T T O 


1 on oTc 


bz.zouz 


6.2 


3.9 


15.92 


1.08 


OTTTTOO 

bU Uoo 


24 


LL2 


1 on noo 

189. 28o 


O CT OO 

DZ.ZOOZ 


8.8 


8.9 


18.94 


2.00 


OT TT TOO 


25 


LL2 


189.284 


62.2539 


37.6 


10.4 


15.83 


4.75 


OTTT TOO 

oU Uoo 


26 


LL2 


1 on noo 

189.288 


DZ.z5z4 


31.9 


17.4 


17.29 


4.67 


OTTT TO O 

bU Uoo 


27 


LL2 


1 on onn 

189. z99 


c o o c oo 


11.0 


5.5 


16.56 


2.58 


OTTTTOO 

bU Ubb 


28 


LL2 


1 on on c 

189.ciU5 


DZ.Z045 


22.8 


8.8 


16.28 


3.33 


O T TT TO O 

bU Uoo 


29 


LL2 


1 on ortr" 

1o9.o0d 


£?0 O/I '70 

Dz.z47z 


32.3 


16.4 


14.91 


3.92 


c T T I roc 
b U LJ ob 


30 


T T O 

LLZ 


1 on Qr\'7 
loy.oU I 


oz.zool 


33.5 


15.4 


17.11 


4.4z 


O T T 1" T O O 

bU Ubb 


31 


LL2 


1 on onn 


i^o o /I no 

Dz.z4Uo 


14.5 


7.2 


19.95 


2.83 


O T TT TO O 

bU Ubb 


32 


LL2 


189.315 


£?0 OOT"? 

02.2377 


14.9 


7.9 


19.77 


2.58 


O T TT TO O 

bU Ubb 


33 


LL2 


ion o nn 

Io9.oz0 


/i^o oono 

oz.ziyz 


15.9 


7.8 


18.03 


3.00 


SUUSS 


34 


LL2 


189.320 


62.2307 


12.1 


9.6 


18.03 


2.83 


SUUSS 


35 


LL2 


189.322 


62.2478 


4.8 


6.2 


15.00 


1.67 


SUUSS 


36 


LL2 


189.324 


62.2328 


6.7 


8.4 


18.12 


1.75 


SUUSS 


37 


LL2 


189.324 


62.2439 


32.5 


12.3 


15.83 


4.75 


SUUSS 


38 


LL2 


189.330 


62.2261 


6.8 


10.2 


18.58 


2.25 


SUUSS 


39 


LL2 


189.333 


62.2269 


24.5 


14.9 


16.65 


2.50 


SUUSS 


40 


LL2 


189.334 


62.2300 


5.4 


4.8 


18.48 


1.17 


SUUSS 


41 


LL2 


189.335 


62.2313 


11.3 


5.7 


15.18 


2.58 


SUUSS 


42 


LL2 


189.340 


62.2291 


24.3 


17.5 


15.37 


2.75 


SUUSS 


43 


LL2 


189.340 


62.2415 


6.9 


5.8 


18.12 


1.67 


SUUSS 


44 


LL2 


189.345 


62.2318 


34.5 


20.3 


15.37 


3.67 


SUUSS 


45 


LL2 


189.346 


62.2385 


21.7 


11.4 


15.09 


3.50 



^ Signal-to-noise ratio estimation when fluxes are integrated over the whole spectral range available 
(14-20^01 for LL2 and 20-35//m for LLl). 

^ Maximum Signal-to-noise ratio achieved in each spectrum and the wavelength at which it peaJjs. 

The spectra were smoothed using a comprehensive banning window over 3 channels in order to 
avoid contribution from any potentially remaining hot pixel value. 

Contracted Spectral bandwidth over which the spectra have flux densities higher than 2a. 

SUUSS 9 is the [SIV] line source discussed in sect. 4.3 which explains the low integrated SNR 
compared to its peak SNR. 
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TABLE 2 

Spitzer Multi- Wavelength Data. 



name 




RAiRs 


DECiRs 


S3.6 " 


54.5 


S5.8 


Ss 


524 ^ 






degrees 


deg 


rees 


fijy 






fiJy 




fiJy 


suuss 


1 


189 


176 


62 


2894 


33.60± 


0.04 


38.60± 0.05 


34.40± 0.30 


23.60± 0.32 


115.0± 5.5 


suuss 


2 


189 


183 


62 


2848 


31.00± 


0.04 


24.70± 0.05 


22.70± 0.30 


41.20± 0.32 


219.0± 5.9 


suuss 


3 


189 


187 


62 


2873 


12.00± 


0.03 


13.40± 0.05 


14.80± 0.30 


lO.lOib 0.32 


162.0± 4.2 


suuss 


4 


189 


191 


62 


2835 


29.30± 


0.04 


22.70± 0.05 


18.20± 0.30 


14.90± 0.32 


80. 6± 4.5 


suuss 


5 


189 


198 


62 


2844 


5.10± 


0.03 


6.35± 0.05 


8.32± 0.30 


6.37± 0.32 


56.4± 4.9 


suuss 


6 


189 


205 


62 


2947 


5.13± 


0.03 


5.66± 0.05 


6.49± 0.33 


6.00± 0.35 


54. 8± 5.4 


suuss 


7 






62 


2764 


18.30± 


0.03 


14.90± 0.05 


ll.SOit 0.30 


20.50± 0.33 


60. 0± 5.9 


suuss 


8 


189 


210 


62 


2827 


10.60± 


0.03 


lO.SOit 0.05 


8.64± 0.31 


7.68± 0.33 


51. 6± 6.8 


suuss 


9 


189 


210 


62 


2869 














f~1 T TT Ton 

suuss 


10 






62 


2936 


3.56± 


0.03 


2.14it 0.05 


0.92it 0.34 


0.54± 0.36 




suuss 


11 


189 


225 


62 


2907 


7.68± 


0.03 


8.02± 0.05 


7.20± 0.34 


6.16± 0.37 


73. 0± 5.8 


suuss 


12 


189 


229 


62 


2829 


45.70± 


0.04 


34.80± 0.06 


28.80± 0.33 


34.10± 0.36 


103.0± 5.7 


suuss 


13 


189 


232 


62 


2822 


16.70± 


0.04 


13.00± 0.05 


11.60± 0.33 


9.09± 0.36 


67.1± 5.0 


suuss 


14 


189 


237 


62 


2788 


6.84± 


0.03 


7.72± 0.05 


9.91± 0.33 


6.69± 0.35 


74.5± 5.7 


suuss 


15 


189 


243 


62 


2638 


12.10± 


0.03 


10.70± 0.05 


7.69± 0.30 


13.90± 0.32 


102.0± 5.5 


suuss 


16 


189 


244 


62 


2634 


12.10± 


0.03 


10.70± 0.05 


7.69± 0.30 


13.90± 0.32 


105.0± 6.1 


suuss 


17 


1 Qn 

loy 




62 


2771 


66.70± 


0.05 


60.20± 0.06 


41.70± 0.34 


120.00± 0.37 


275.0± 5.7 


suuss 


18 


189 


246 


62 


2610 


13.20± 


0.03 


13.20± 0.05 


11.60± 0.29 


7.69± 0.32 


47.2± 6.3 


suuss 


19 


189 


251 


62 


2687 


15.70± 


0.03 


11.30± 0.05 


8.33± 0.32 


8.69± 0.34 


103.0± 8.5 


suuss 


20 


1 Qn 


Zol 


62 


2714 


36.60± 


0.04 


27.20± 0.06 


29.90± 0.33 


26.50± 0.35 


502. 0± 7.7 


suuss 


21 


189 


259 


62 


2532 


15.50± 


0.03 


10.60± 0.05 


8.76± 0.29 


7.55± 0.31 


71. 0± 6.2 


suuss 


22 


189 


263 


62 


2519 


6.58± 


0.03 


5.94± 0.05 


2.66± 0.29 


3.52± 0.31 




suuss 


23 


189 


275 


62 


2552 


13.20± 


0.03 


9.33± 0.05 


9.37± 0.31 


7.66± 0.33 


59. 2± 6.3 


suuss 


24 


189 


283 


62 


2582 


8.79± 


0.03 


11.50± 0.05 


14.30± 0.32 


10.70± 0.36 


219.0± 6.5 


suuss 


25 


189 


284 


62 


2539 


27.40± 


0.04 


18.80± 0.05 


20.10± 0.32 


14.80± 0.35 


198.0± 5.9 


suuss 


26 


189 


288 


62 


2524 


26.40± 


0.04 


24.50± 0.06 


18.50± 0.32 


17.90± 0.35 


123.0± 5.4 


suuss 


27 


189 


299 


62 


2538 


21.50± 


0.04 


18.70± 0.06 


13.30± 0.34 


29.90± 0.38 


88. 1± 6.2 


suuss 


28 


189 


305 


62 


2548 


16.50± 


0.04 


12.00± 0.06 


12.30± 0.35 


9.93± 0.41 


169.0± 6.3 


suuss 


29 


189 


306 


62 


2472 


74.10± 


0.05 


54.80± 0.06 


46.30± 0.34 


33.00± 0.38 


241.0± 7.0 


suuss 


30 


189 


307 


62 


2531 


34.70± 


0.04 


27.70± 0.06 


26.30± 0.35 


32.80± 0.41 


149.0± 6.2 


suuss 


31 


189 


309 


62 


2403 


31.50± 


0.04 


25.40± 0.06 


21.70± 0.33 


21.80± 0.36 


103.0± 6.2 


suuss 


32 


189 


315 


62 


2377 


16.30± 


0.04 


ll.OOit 0.06 


12.20± 0.33 


7.94± 0.37 


76. 6± 6.5 


suuss 


33 


189 


320 


62 


2292 


2.32± 


0.03 


1.68± 0.05 


1.34± 0.31 


0.66± 0.35 




suuss 


34 


189 


320 


62 


2307 


18.00± 


0.04 


15.60± 0.05 


ll.lOit 0.32 


12.60± 0.36 


66. 0± 5.3 


suuss 


35 


189 


322 


62 


2478 


10.20± 


0.04 


7.19± 0.06 


7.06± 0.36 


5.71± 0.42 


81. 1± 6.4 


suuss 


36 


189 


324 


62 


2328 


14.20± 


0.04 


11.90± 0.06 


8.11± 0.33 


16.70± 0.37 


54. 5± 5.6 


suuss 


37 


189 


324 


62 


2439 


22.50± 


0.04 


20.80± 0.06 


15.50± 0.36 


38.50± 0.42 


119.0± 6.7 


suuss 


38 


189 


330 


62 


2261 


1.33± 


0.03 


1.38± 0.05 


0.833± 0.31 


1.87± 0.36 




suuss 


39 


189 


333 


62 


2269 


77.90± 


0.05 


61.30± 0.06 


50.60± 0.33 


63.80± 0.37 


185.0± 4.6 


suuss 


40 


189 


334 


62 


2300 


35.50± 


0.04 


25.50± 0.06 


23.00± 0.34 


17.30± 0.39 




suuss 


41 


189 


335 


62 


2313 


35.50± 


0.04 


25.50± 0.06 


23.00± 0.34 


17.30± 0.39 


142.0± 5.6 


suuss 


42 


189 


340 


62 


2291 


43.50± 


0.05 


34.90± 0.06 


26.70± 0.34 


25.60± 0.39 


217.0± 6.4 


suuss 


43 


189 


340 


62 


2415 


6.52± 


0.04 


5.31± 0.06 


4.34± 0.38 


3.01± 0.43 




suuss 


44 


189 


345 


62 


2318 


2.60± 


0.04 


1.79± 0.06 


1.17± 0.35 


1.66± 0.42 


301.0± 4.7 


suuss 


45 


189 


346 


62 


2385 


32.60± 


0.04 


24.90± 0.07 


19.20± 0.38 


16.80± 0.43 


148.0± 5.8 



^ IRAC photometry from GOODS legacy progra m, jPickinson et aT] ( |2003^ 
^ MIPS 24um photometry from GOODS-North, ICharyl ||2007|l 
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TABLE 3 
Templates IR properties. 



Template 


z 


6.2/ini EQW 
fim 


7.7^"! EQW 
fj,m 


log{LiR[LQ]) 


Ref. 


Mrk^ol 


0.042 


0.011 




12.57 


1 


TV T 1 o ^ O 


0.038 


0.171 




12.15 


1 


Mrk4Do 


0.051 




0.029 


11.70 


2 


TV 1 n 1 /I 

Mrkl(J14 


0.163 


0.048 


0.135 


12.53 


2, 7 


JNCjOdz4U 


0.025 


0.52 


2.60 


11.85 


3 


T T/^ /~i r 1 r\-i 


0.039 


0.188 


0.419 


12.00 


1, 2 


Arp220 


0.018 


0.253 




12.16 


1 


TT~> A o n 1 orv 


0.042 


0.035 




12.16 


1 


IRAS 08572 


0.058 


<0.012 




12.14 


1 


TT3AQ 101T0 


0.073 


0.517 


0.569 


12.33 


1 


IRAS 14348 


0.083 


0.254 




12.35 


1 


IRAS 15250 


0.055 


0.023 




12.05 


1 


IRAS 22491 


0.077 


0.594 


0.671 


12.19 


1 


NGC1569 


~0 




0.402 


8.76 


4 


3C120 


0.033 








5 


PG1612+261 


0.131 








5 


PAH 1^ 










6 


PAH 2=^ 










6 


PAH 3=^ 










6 


PAH 4=^ 










6 


PAH 5=^ 










6 



References. — ( 1) lArmus et al l j2007|l; (2) |Armus et al] l |2004[ |; (3) 
|Armus et al.|||2 006|; f4)|Wuef^r^ll^nD B|; (5) K Ogle , private communication; 
(6) [Smith efa l. (2007b); (7) Kim & S andersj ( [l998t . 

As explained in Smit h et al.| ( |2007b^ , those 5 templates are composite spec- 
tra computed in arbitrary umts ot ul^. We thus do not have corresponding 
redshift, luminosity or equivalent width. 
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TABLE 4 
Cross-correlation analysis results. 



ID 




Field 






2re/ type 


template 


spectral type ^ 


QTTTTQQ 


1 


I 


1.70 


1.59 


p 






suuss 


2 


1 


0.50 


0.50 


S 


PAH4 


PAH 


CT TT TQQ 


3 


1 


2.03 


2.03 


b 


Arp220 


mixed 


QT TT TG Q 


A 

4 


i 


u.yo 


i.Ul 


r 


TD A U/l 


T> A u 


QT TT TGC 


Q 
O 


i 


n fit; 






TD A TTi-> 


Tr> A XT 


CT TT TQQ 


9 


1 


2.08 


■ ■ ■ 




iNLjOlooy 


line 


CT TT TCQ 
bU Ubb 


14 


1 


2.20 


2.67 




TV (T».KO'7Q 

Mrkz / v> 


mixed 


CT TT TGC 
bU Ubb 


1 


i 


U.41 


U.4u 


b 


Jr^Arlo 


Tr> A u 


CT TT TCC 
bU Ubb 


17 


1 


0.30 


0.30 


b 


pglolZ+ZDl 


line 


CT TT TCC 
bU Ubb 


18 


1 


2.04 






TT3AC ic;oc;n 


blU 


CT TT TCC 
bU Ubb 


on 


i 


U.yl 


n ni 

u.yi 


b 


Tr> A U/l 


T> A tJ 
Jr^Arl 


CT TT TO C 

bU Ubb 


24 


2 


2.03 


2.15 


p 


PAH5 


PAH 


CT TT TC C 

bU Ubb 


25 


2 


0.84 


0.84 


s 


PAH4 


PAH 


CT TT TC C 

bU Ubb 


26 


2 


1.20 
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^ Spectroscopic (S) redshifts are from the Treasure]) Keck Redshift Survey (Wirth et 
al., 2004) and photometric (P) redshifts are from Caputi et al. (2006). 
^ Sort the main spectral signature identified in the MIR spectra: PAH for strong 
aromatic features, SiO for silicate absorpt ion, mixed for intermediate spectra, and 
line for the 2 line sources discussed in Sect 14.31 

Sources with no conclusive identification from the cross-correlation method but for 
which we had ancilliary redshifts (same reference). 
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